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SUM&RY 


REVIEW OF PREVIOUS EXPERIMENTAL WORK 


The drift field model used by Wolf to calculate 
the short-circuit current (l sc ) was extended to permit 
calculations of the open-circuit voltage (V Q _) and the 
maximum power (P.) under conditions of lUunt nation of 
either tungsten (2300° C) source or AMD sunlight. 
Voltages were calculated using an expression for the 
drift field diode saturation current derived here. 

The model, applied to the oxygen rich (C -13 group) 
lithium solar cells, was used to calculate the pre- 
and poet-electron bombardment trends of the V QC , P m , 
and I 8C for lithium gradients in the range of 
10*" to 10*9 Li/cnr. Published experimental data 
characterizing these cells were used to tailor the 
model. The calculated trends are in reasonable agree- 
ment with the empirical data of Faith. Diffusion 
length degradation and carrier removal effects were 
sufficient to predict the cell performance up to 
3 x 10 ** e/anr. Beyond this fluence it was necessary to 
include drift field removal effects. 

INTRODUCTION 

The unusual performance of lithium- containing 
solar cells has been studied by several investigators 
(1-3). In order to explain the self-healing effects 
of these cells when exposed to a radiation et /ironment , 
it has been proposed that the lithium neutralizes the 
radiation induced defects resulting in less severe 
effect on cell performance. Recently Faith (4, 5) has 
measured , using a tungsten light source, the open cir- 
cuit voltage ( V oc ) , short circuit current ( l S r ) , and 
maximum power (P m ) of oxygen rich (C-13 group) lithium 
cells exposed to 1 MeV electrons in the 3X10*3 to 
3Xl<r*5 e/cnr fluence range. From these data, empiri- 
cal formulae were constructed which relate the radia- 
tion behavior to the magnitude of the lithium gradient 
near the junction. Additional work by lies (6) had 
shown that unbombarded cell performance can be corre- 
lated with the lithium gradient. In addition Faith 
et al. ( 7 ) have considered a case in which the lithium 
gradient gives rise to an Internal electric field. 

Although a drift field model has been developed 
and applied to calculate the short circuit current 
performance ( 8 - 10 ), the model has not been used to 
analyze the open circuit voltage or the maximum power 
of the lithium cell. Thus it is the purpose of this 
paper to formulate a general field model for the lith- 
ium solar cell, to use this model to calculate the 
current, voltage, and power of the cell, and to pre- 
dict the initial and radiation bombarded performance. 
The model will be applied specifically to the C-13 
group of lithium cells, since the electrical perfor- 
mance of these cells is well documented. The formula- 
tion for Ig C obtained by Wolf ( 8 ) and the diode sat- 
uration current derived here will be used to obtain 
the V oc from the simple diode equation. The maximum 
power and cell curve factor are also calculated. 

These equations include drift field effects and apply 
to the low injection case. 


Lithium Gradient ana Cell Performance 

The calculations using the four layer drift field 
model are based on measurements of the lithium profile 
and diffusion length for the C-13 group of oxygen-rich, 
lit hi\m-d if fused p an n solar cells published by Faith 
(4, 5 ) and lies ( 6 ). These data are summarized in 
Thble I for nominal lithium concentrations evaluate 
at the junction edge (x* ) , at the rear contact (W) , 
and at the subregion edge (X). The pre irradiated cell- 
performance, e.g. , V F ffi , I gc , and lithium rofile, 
was measured for cells In this group and was found to 
correlate with the lithium gradient at the junction 
edge. The high vpltage cells, wnicn nad a high gra- 
dient (~ 10*9 cm**), were fabricated in short lithium 
diffusion times (375° C - 180 min). Using experimental 
data, it was estimated that the gradient extended over 
a region 3 uffi thick next to the junction. Similarly, 
cells which were diffused for longer times (325° C - 
480 min) had a lower gradient (~ 10 *" an*) and also a 
lower Vjxj compared to the high gradient cells. The 
subregion thickness for this case was estimated to be 
about 10 to 20 urn. Cells diffused in the short times 
had higher measured lithium concentrations everywhere 
in the n-base, than were cells diffused for longer 
times. No experimental diffusion schedule can be 
as. igned to the medium gradient ( 5 <10*" ) cells. How- 
ever, it is thought that these gradients are a result 
of the 373 ° - 180 min diffusion and extend over a 
thickness that is intermediate to the high and low 
gradient cases. Correlations of tte. lithium gradient 
with the measured values of I gc and were also 
made. In these cases however the trend appears to be 
that the highest values of I gc and P occur for the 
low gradient and the lowest values for the high gra- 
dient. The carrier lifetimes shown in Table I repre- 
sent nominal data representative of each type cell. 

Empirical Bombardment Data 


Measurements of the V QC , P , and I of select- 
ed C-13 cells after exposure to 1 MeV electron fluences 
in the range of 3*10*3 to 3*10*5 e/cm 2 have been re- 
ported by Faith (4,5). A tungsten light source was 
used for cell evaluation. From these data empirical 
formulae were developed that reflect a dependence of 
each cell parameter on the electron fluence and on the 
lithium gradient (dN/dx). These relationships which 
are valid only in the above fluence range are shown 
below: 


Open Circuit Voltage 

V 0C = v(o)| $ + A log (10-tf g) 

Maximum Power 

P m . P( 0 )| 4 - B log (lO-“ g ) 


( 1 ) 

( 2 ) 


l 





Short Circuit Current 


P 


( 6 ) 


i 8C - i(o)| 4 - c log (wr»g) (3) 

The slopes A, B, and C, determined from the measured 
data, are 16 mv/decade , 2.4 mw/decade, and 8.5 na/ 
decade respectively. The intercepts V(o)|^ , P(o)L 
and J(o)U represent values corresponding to the 9 
lxiolS cm* 1 * gradient and have different magnitudes at 
eacn fluence level. In addition an empirical relation- 
ship for the 1 MeV electron damage coefficient (K L ) 
was determined from the I fiC bombardment behavior and 
is shown below. 

•4,(0) = (^) (1 - 0.063 log *) (U) 


CF X V XI 
oc s 


The curve factor (CF) v&a calculated using the 
equations and an iteration method proposed by Infers hi 
(11). The overall computer program is arranged to 
permit performance analysis of a variety of cell de- 
signs and operating environments. The cell geometry, 
the depletion region width assuming an abrupt junction, 
material properties and drift field strength can be 
selected or calculated in addition to the choice of 
cell operating temperature and photon spectrum for 
either tungsten or sunlight. The cell spectral re- 
sponse can also be calculated. Reflection coefficients 
for either a bore or SiO coated cell can be selected 
to simulate the real cell, rather than assuming com- 
plete absorption. 

Modeling of he Lithium Cell 


THEORY 

Description of 4- layer Model 

The cell model shown in figure 1 is identical to 
the one used by Wolf (8) to calculate drift field cell 
collection efficiencies. The homojunctlon semiconduc- 
tor device Is divided into a diffused and a base re- 
gion. Each region is further divided into two sub- 
regions resulting in four layers for the total struc- 
ture. The subregion boundaries, Y and X, can be varied 
to alter the relative thickness of each subregion. The 
depletion region is not counted as a fifth layer since 
all photogenerated carriers in this region are assumed 
to be collected (i.e. , 100% collection efficiency). 

The impurity concentrations shown In figure 2 are 
specified at the boundaries of the front and rear 
surface (0 and W ) , the subregions X, Y and the deple- 
tion region. This arrangement permits drift electric 
fields to be applied across the subregions. For this 
study only exponential impurity distributions are 
assumed which result In constant drift fields. Other 
cell properties used in calculations are shown in 
Table II. Individual subregion properties such as 
mobility (u) and carrier lifetime (t) can be specified 
independently. Additionally the surface recombination 
velocity for the rear (s r ) and front (s^ ) surfaces are 
specified. The calculations are for a tungsten 2800 K 
spectrum except where noted. 

The general approach to the solution of this 
model requires solving the continuity equation using 
the current transport equation and appropriate boun- 
dary conditions. 

We have extended the work of Wolf (8) to permit 
evaluation of the open circuit voltage and maximum 
power, in addition to the short circuit current. Using 
the seme boundary conditions at the subregion and front 
and rear surfaces, as cited in (8) an expression for 
the drift field diode saturation current (Iq) was de- 
rived, and is shown in Appendix A. Hence Iq reflects 
the four layer concept in that not only can the mate- 
rial properties of each sublayer be specified, but also 
that drift fields can be introduced via the exponen- 
tial inpurity distributions. Thus expressions for Iq 
and I 8C tere programmed on a high speed digital com- 
puter and used to calculate the voltage from the ex- 
pression 



The maximum power was calculated using the well known 
relationship 


The computational model used to represent the 
experimental cells was constructed in the following 
way. The lithium concentration from figure 2 is spe- 
cified at the depletion region edge (x * *.), at the 
subregion boundary, and at the rear contact X e W. 

Even though the real profile may be .far more complex, 
the lithium concentration was approximated by an ex- 
ponential increase from the depletion edge to the sub- 
region boundary. From there to the ohmic rear con- 
tact , the lithium was assumed uniformly distributed 
and therefore field free. Thus, having specified the 
concentrations, the width of the subregion nearest 
the junction (l.e., field region width) was calculated 
using the experimental values for the lithium gradient 
and equation (B6) derived in Appendix B. (Conversely, 
given the concentrations and the field region width, 
the gradient could be calcula ed if so desired. ) The 
field region widths calculated for the l<10 i 9 and 
lxio^ cm* 1 * lithium gradients are about 3 and 11 urn 
respectively. These values are in reasonable agree- 
ment with the experimental c'au. 

The selection of values for the various cell 
parameters also was based, where possible, on experi- 
mental data. The initial carrier lifetimes shown in 
Table I were based on the I 8C versus diffusion length 
data from Faith (12). The same lifetime was assumed 
for both subregions. The differing carrier mobilities 
shown in Table II reflect the difference in average 
lithium concentration of the two subregions. Other 
cell specifications used in the modeling are sham in 
Table II. 

For the electron fluence range of zer:> to 3xl0 11+ 
e/cn£, the lithium cell performance was easily simu- 
lated assuming only diffusion length (lifetime) degra- 
dation and carrier removal effects. The changes in 
the initial diffusion length were calculated using the 
damage coefficients computed from equation (4), Car- 
rier removal effects were simulated by reducing the 
initial lithium concentrations using a constant carrier 
removal rate (CRR) of 0.04 carriers/cm- electron. 

At the 3x10^5 e/cm 2 fluence the assumption of 
field removal, in addition to diffusion length and 
carrier removal effects, was necessary to simulate the 
measured data. Field removal was included by artifi- 
cially reducing the electric field strength in the 
subregion nearest the junction. The field was de- 
creased using a trial and error method rncil the vol- 
tages calculated for the 3X10*5 e/cm 2 fluence were 
brought into agreement with the experimental data. It 
is obvious that this procedure "forced” the highest 
fluence voltages into agreement with the experimental 
trends. However the appropriateness of this approach 
was judged primarily by the magnitudes and trends that 
resulted for the short-circuit current and the maximum 
power. 
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Open- Circuit Voltage 


RESULTS 


The results of the calculations for the Initial 
and bombarded v oc» p max and I gc are shown in 
figure 3. fQr Hthl\m gradients In the range of 10" to 
10*9 cm . The calculated data are based on lifetimes 
calculated from the empirical damage coefficients and 
for a constant CRR of about 0.0k. Figure 3a shows the 
calculated voltage data as the open symbols compared 
to the empirical data, shown by the solid lines. The 
initial calculated voltages are representative of the 
range (error bars) of experimentally measured values 
corresponding to gradients of 10*8, 5*lCr^, and 10*9 
an"^. In general it was found that the initial vol- 
tages depend strongly on the lithium concentration at 
the depletion region edge, and to a lesser degree on 
the initial carrier lifetime. The agreement between 
the calculated and experimental trends at the three 
fluence levels is excellent. The fall off of voltage 
being about l£ mV/ decade; in agreement with experiment . 

Decreases in the lithium gradient are predicted 
above 3 X 10*^ e/cm and are shown in the figure. TLese 
changes are due to carrier removal effects and are 
most severe for the lowest gradient due to an overall 
lower lithium concentration in that cell. It was 
found that lifetime degradation and carrier removal 
alone would result in a slope at the 3x10*5 e/cm 2 
level that was much steeper than that indicated by 
the expe Imental data. Agreement was achieved, how- 
ever, by artificially reducing the drift field to 
simulate radiation induced field removal. Although 
there is no experimental evidence suggesting that 
field removal does indeed occur, Dresselhaus (13) sug- 
gests that mobile negatively charged radiation induced 
defects (1, 14) can drift within the field region. 

If this were the case a counter field would arise 
weakening the lithium drift field. 

Maximum Power 

Figure 3b shows the initial and bombarded trends 
predicted for the maximum power. The calculated curve 
factors of about 0.80 were reduced by 8$ to more 
closely agree with experimental data. These data 
correspond to the same conditions as discussed for the 
voltage. The solid lines represent the empirical data, 
while the broken lines and open symbols show the cal- 
culated trends. The experimental data were converted 
to a power density by assuming an active area of 1.8 
cm 2 . The disagreement between the experimental and 
the calculated power densities is about *5 to 10$. 

The predicted slopes are slightly steeper but yet In 
reasonable agreement with the empirical trends. As 
shown in Table III, the calculated range is 3*33 to 
2.4 6 mW/decade compared to the average experimental 
value of 2.40 mW/decade. 

Short-Circuit Current 

The predicted short-circuit current densities 
(broken lines and open symbols) are shown in figure 3c > 
and are again compared with the empirical trends 
(solid lines). The experimental current densities 
were calculated assuming an active area of 1.8 os 2 . 

As can be seen the predicted dependence of current on 
lithium gradient is in reasonable agreemant with ex- 
periment. These slopes are compared in Table III. 

The calculated values range from 6. 91 to 8.37 ®a/ 
decade compared to the experimental average of 8.5 
ma/decade. It is perplexing, however, that the magni- 
tudes of the calculated current density at the 3x10*3 
and 3x10*5 e/cm 2 fluencea are in disagreement by about 
10$ with the measured data. The current is under- 


estimated at 3X10*3 e/c m 2 , but overestimated at the 
3x10*5 fluence level. 

Since the current is more sensitive than the V QC 
or P ffl to intensity, spectrum, and diffusion length, 
failure to closely simulate the experimental conditions 
would have a more pronounced effect on the current. 
Furthermore, the radiation damage trends are found to 
be sensitive to the character cf the light source. 
Figure 4 shows the results of the I gc calculation 
where the AM) sunlight spectrum was substituted for 
the tungsten 2600° K spectrum used previously. Not 
only is the slope significantly shallower ( 4.8 ma/ 
decade) than for the tungsten case, but the AMD current 
densities are also about l£$ higher. These results are 
not unexpected since it is well known that a radiation 
damaged cell will appear more severely degraded when 
evaluated with a red-rich light source compered to a 
similar evaluation with a blue peaked source. This 
8 hows that in addition to the intensity, the spoctrur 
of the light source plays a significant role i:. deter- 
mining a cell's radiation performance trends as a func- 
tion of lithlimi gradient. 

CONCLUSION 

The four layer drift field model was used to pre- 
dict the pre- and pest -irradiation performance of the 
lithixm solar cell. Using this model and published 
experimental data for the C-13 group of lithium cells 
the V oc , I gc , and P m were calculated and compared 
with the measured performance. The empirically deter- 
mined dependence of these parameters on electron 
fluence and lithium gradient has been verified using 
this model and by assigning diffusion length degrada- 
tion, and carrier removal. For flue nee e above 3*10*^ 
e/c nr the added assumption of field removal was made. 
Calculations using an AMO light source epectrum pre- 
dict that th.. rate of change of the performance para- 
meters with lithium gradient is less than when measured 
under tungsten illumination. 

It is concluded that this model can be applied 
generally ^ explain the performance of all lithium 
cells providing sufficient information is available 
concerning the lithium profile and the minority carrier 
lifetime, as well as the character of the light source. 

APPENDIX A - SATURATION CURRENT 


Using the model shown in figure 1, the n-base 
component of the diode saturation current density is 
derived from the continuity and the current transport 
equations. The treatment Js identical for the p-base 
component. Hence, for the one dimensional case, the 
general differential equation can be written in the 
steady state 


in qE 



(Al) 


where £ is the constant drift field and 1^ the 
n-base bulk diffusion length. Equation (Al) has the 
general solution 

n(x) « D x e riX + Dg e r 2 x (A2) 

where 
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rp is similar except that the 2nd term has a negative 
sign instead of positive. A similar solution with 
different coefficients is written for the second sub- 
region. Squat ion (A2) is the general expression for 
the minority carrier density in the n-base. The drift 
field solution is Obtained by solving for the coeffi- 
cients subject to the following boundary conditions: 


A similar eqmtion can be derived for the diffused 
layer. Therefore the total diode saturation current 
used to calculate tne open-circuit voltage can be 
written 

l o m \*\ (A5) 


1. at x s 

4 

n(xj) = ~ (exp 3V-1) 

where 9 = q/kT, n^ the intrinsic concentration, and 
Np the lithium concentration at the junction edge. 

2. at x = X 

rijOO = a,(X) 

dn,(X) diu(X) 

D n, S~ + = D + uK^X) 

(A3) 

where the subscripts 1 and 2 refer to the subregion 
nearest the junction and nearest the contact respec- 
tively. 


APPENDIX B - FIELD REGION WIDTH 

The field region width was calculated by assuming 
lithium concentration (N2 ) at the depletion edge (xj) 
tlat increases exponentially to the subregion boundary 
X. Tte width of this subregion is therefore (X-x*). 

Tte concentration (%) is assessed constant from tne 
subregion boundary (x » X) to the rear contact (x * W). 
The lithium concentration within the subregion 
(X - x j ) can be expressed 

N^(x) « exp ax (Bl) 

where a is the grade constant. (Bl) is subject to 
the boundary condition: 

at x * Xj: N( x j) * N o ex P * n l> C 02 ) 

at x = X: N(X) = N exp aX * N 
o n 


3. at X S W 

dn^w) 

dx 


- S^V) 


where 


_ f2 Jr 

S 2 * KT + D_ 


and s r is the surface recombination velocity. 

After an algebraic manipulation, the following 
expression is obtained for the diode saturation 
current 

\ '"“‘A r 


where 


c+ espfel- A l)s 
1/2 




•!&-[(%) ♦(£)*) 


s F -5 

nm nm (AU) 


F + G 
nm nm 


where m « 1 is the subregion nearest the junction, 
m * 2, nearest the contact ; and 


(- 3^}*P 20n2(H ‘ x) [>V A l] * fv*i] 




The lithium gradient is given by 

~ m aN exp ax 
dx o ^ 

evaluating (B3) at the depletion region edge x c x^ 


(B3) 


dN 

dx 


a n t 


(BU) 


From (B2), tte grade constant can be determined 


In 


X- x. 


(B5) 


Substituting (B5) into (B4) and solving for X - x^ 


X- Xj -N i: 


N, 

H 


(B6) 
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High 

Gradient 


Medium 

Gradient 


Low 
Gradient I 


Field Region (X-x-j ) 
Junction Depth (xj) 

Cell Thickness 

(ym) 


3.0 


6.6 
0 . 25 * 
• 375 < 


11.1 


Mobility, 

jcm^/volt-sec 


(X-Xj) Region 
(W-X) Region 


300 


. 500 ■ 
300 


400 


Diffusion Region Profile 


erfc (1/3 Dead layer) 


Surface 

| Recombination] 
Velocity, 
cm/sec 


Front 

Rear 


' 10 5 


■ 10“ 


AR Coating 


• SiO ■ 


Lit him Gradient, 
atoms/cnr 


l.lxiO 1 ? 




10 18 


9 . 5 - 1 C 17 I 


11. J. J. Loferski, "Theoretical Considerations 
Governing the Choice of the Optimum Semiconductor 
for Photovoltaic Solar Energy Conversion,” 

Journal of Applied Physics, Vol. 27 , No. 7 , 

July 1956. 

12. T. J. Faith, Jr. , "Density and Fluence Dependence 
of Lit him Cell Damage and Recovery Characteris- 
tics," JPL, April 29 , 1971 , p. 13 . 

13. M. S. Dresselhaus, "Radiation Damage Annealing 
Kinetics," ibid, p. 101-113. 

14. G. 0 . Watkins, J. Phys. Soc. Japan 18, 22 ( 1963 ). 


TABIE I 

PUBLISHED EXIER RENTAL DATA FOR C-13 GROUP 
OF LITHIUM CELIE - REFS. 5 AND 12 











x=x j 

8> ao 11 * 

3 X 10 11 * 

Lithium 


Concentration^ 

x=X 

5KL0 16 

1 X 10 16 

atoms/ cm3 


5 x 10 ^ 

1*10^ 

Gradient, atoms/czn* 

IxkP-9 

lxlO 1 ^ 

Initial Lifetime 
(N-base)( 2 ', usee 

10 

27 

Subregion Width(3) >Mm 

-3 

10-20 


TABIE III 


CALCUIATED AND EXIERRCNTAL RADIATION 
DAMAGE DATA FOR Ig CC AND FOR 

HIGH AND LOW GRADIENT LITHIUM 
DRIFTED SOIAR CELTS 


Current (Igcc) 

Electron 

Al (Ael.8 cm 2 ) 

Experimental 

e/cm 2 

Decade 

Slope* 

3 XI 0 13 

8.46 

8.4 


8.91 

8.4 

3 X 10 15 

8.37 

8.9 

Average 

8.4 

8.5 

Maximum Power 

Electron 
Fluence , 
e/cm 2 

AP (Aal.8 an 2 ) 

Experimental 

Decade 

Slope* 

3*10 13 

3.33 

2.39 


3 . 21 * 

2.40 

3*L0 15 

2.46 

2.1*2 

Average 

3.01 

2 . 1*0 


* 


Refs. 4 and 5 
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CROSS SECTION OF FOUR LAYER MODEL 
OF LITHIUM SOLAR CELL 
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SUBREGION VARIABLES: E, T, p 


SUBREGION BOUNDARIES 
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METAL CONTACT 




Fig. 1 


LITHIUM PROFILE AND FIELD REGION WIDTHS 


DONOR 

CONCENTRATION 

(CM' 3 ) 
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OPEN CIRCUIT VOLTAGE IMMEDIATELY AFTER IRRADIATION 
VERSUS LITHIUM GRADIENT FOR TUNGSTEN 2800K SPECTRUM 



i 
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MAXIMUM POWER IMMEDIATELY AFTER IRRADIATION AGAINST 
LITHIUM GRADIENT FOR TUNGSTEN 2800K SPECTRUM 



Fig. 3b 
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SHORT CIRCUIT CURRENT IMMEDIATELY AFTER 
IRRADIATION AGAINST LITHIUM GRADIENT 
FOR TUNGSTEN 2800K SPECTRUM 



CALCULATED RADIATION BEHAVIOR OF SHORT CIRCUIT 
CURRENT DENSITY FOR AIR MASS ZERO SPECTRUM 



Fig. 4 
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